Post-Newtonian constraints on f(R) cosmologies in metric formalism 
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We compute the complete post-Newtonian limit of the metric form of f(R) gravities using a 
scalar-tensor representation. By comparing the predictions of these theories with laboratory and 
solar system experiments, we find a set of inequalities that any lagrangian f(R) must satisfy. The 
constraints imposed by those inequalities allow us to find explicit bounds to the possible nonlinear 
terms of the lagrangian. We conclude that the lagrangian f(R) must be almost linear in R and that 
corrections that grow at low curvatures are incompatible with observations. This result shows that 
modifications of gravity at very low cosmic densities cannot be responsible for the observed cosmic 
speed-up. 
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I. INTRODUCTION 

It is now generally accepted that the universe is 
undergoing a period of accelerated expansion pj 0, 
which cannot be justified by the description provided 
by the equations of motion of General Relativity (GR) 
and a universe filled with standard sources of matter 
and energy. It has been suggested that this effect 
could have its origin in, among other possibilities, 
corrections to the equations of motion of GR generated 
by nonlinear contributions of the scalar curvature in 
the gravity lagrangian f(R) @, 0, 01 (see also @). 
Reasons for considering nonlinear curvature terms in 
the gravity lagrangian can be found in quantum effects 
in curved space Ml or in certain low-energy limits of 
string/M-theories |8j ■ The nonlinearity of the lagrangian 
can also be related to the existence of scalar degrees 
of freedom in the gravitational interaction 0. In any 
case, the fact that certain f(R) lagrangians naturally 
lead to early-time inflationary behaviors is the main 
motivation to study possible new gravitational effects in 
the late-time cosmic expansion. 

Once a nonlinear lagrangian f(R) has been proposed, 
the equations of motion for the metric can be derived 
in two inequivalent ways. On the one hand, one can 
follow the standard metric formalism, in which vari- 
ation of the action with respect to the metric leads 
to a system of fourth-order equations. On the other 
hand, one may assume that metric and connection 
are independent fields and then take variations of the 
action with respect to the metric and with respect to 
the connection. In this case, the resulting equations 
of motion for the metric are second-order. Only when 
the function f(R) is linear in R, GR and GR plus 
cosmological constant, metric and Palatini formalisms 
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lead to the same equations of motion. In this work we 
will be exclusively concerned with the metric formalism. 
The Palatini formalism will we considered elsewhere ^3 • 

Though much work has been carried out in the last few 
years with regard to f(R) gravities in the cosmological 
regime, very little is known about the form that the 
gravity lagrangian should have in order to be compatible 
with the most recent cosmological observations [TT| . 
The main reason for this seems to be the fact that the 
precision of the supernovae luminosity distance data and 
other currently available tests supporting the late-time 
cosmic speed-up is not enough to discriminate with 
confidence between one model or another. It would be 
thus desirable to have a new arena where to test these 
theories with higher precision. In our opinion, the solar 
system represents a scenario more suitable than the 
cosmological one to study the possible constraints on 
the lagrangian f(R). In fact, if in addition to modified 
gravitational dynamics, sources of dark energy were 
acting in the cosmic expansion, it would be very difficult 
to distinguish their effect from a purely gravitational 
one. In the solar system, however, it is ordinary matter 
which dominates the gravitational dynamics, being the 
contribution of dark sources negligible. Therefore, we 
should see the solar system as a suitable laboratory to 
impose the first useful constraints on f(R) cosmologies. 

In order to confront the predictions of a given grav- 
ity theory with experiment in the solar system, it is 
necessary to compute its weak-field, slow-motion (or 
post-Newtonian) limit. This limit has been computed 
for many metric theories of gravity and put in a stan- 
dardized form |12| , which depends on a set of parameters 
that change from theory to theory (Parametrized 
Post-Newtonian [PPN] formalism). However, for f(R) 
gravities this limit has not yet been computed in detail. 
The Newtonian limit of these theories was recently 
studied in |rj ■ On the other hand, since the metric form 
of f{R) gravities can be represented as the case to = 
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of Brans-Dicke-like scalar-tensor theories, it is tempting 
to use the post-Newtonian limit of those theories given 
in the literature 0, Q] (see also 0|) to check the 
viability of particular models. This was proposed in 
fill , where it was concluded that the Carroll et al. 
model Q (in metric formalism) was ruled out according 
to the observational constraints on the parameter 7 
corresponding to that model. That result was based on 
the fact that the scalar field had a small effective mass, 
which was computed in terms of the second derivative 
of the potential. However, that prescription is usually 
derived under the assumption that the potential and 
its first derivative vanish (sec for instance [Til 
conditions that, in general, cannot be imposed on f(R) 
theories (see section ITTTf) . 

In this work we will use the scalar-tensor representa- 
tion introduced in |l6j to compute the post-Newtonian 
limit of f{R) gravities taking into account all the terms 
associated to the potential of the scalar field. In other 
words, we will not make any assumption or simplification 
about the function f(R) that defines the lagrangian. 
We will actually compute the post-Newtonian limit 
corresponding to Brans-Dicke-like scalar-tensor theories 
with arbitrary potential and a generic constant value of 
lj and will then particularize to the case u> = 0, which 
corresponds to the metric form of f(R) gravities. In this 
manner we generalize the results of the literature so as to 
include all the terms that are relevant for our discussion. 
The Palatini form of f(R) gravities, which can be 
represented as the case oj — —3/2 of Brans-Dicke-like 
theories, represents an exception of the general case lo = 
constant and will be studied elsewhere Hoi. 



The resulting post-Newtonian metric will allow us 
to confront the predictions of these theories with the 
observational data. In this way, we will find a series 
of constraints for the lagrangian, which is a priori 
completely unknown. Those constraints turn out to 
be so strong that the lagrangians compatible with 
observations are bounded by a well defined function 
that prohibits the growing of the nonlinear terms at low 
curvatures. This result will be enough to invalidate the 
arguments supporting the cosmic speed-up as due to 
new gravitational effects at low curvatures. 

The paper is organized as follows. We first derive 
the equations of motion in the original f(R) form and 
show how to obtain the scalar-tensor representation out 
of them. Then we comment on the choice of coordinates 
and boundary conditions. The post-Newtonian limit is 
computed in the Appendix and commented in section llVl 
where we obtain constraints on the lagrangian from the 
experimental data. Section[V]is devoted to the discussion 
of particular models. In section IVII we find the form of 
the lagrangian that satisfies the constraints. We conclude 
the paper with a brief summary and conclusions. 



II. EQUATIONS OF MOTION 

The action that defines f(R) gravities in the metric 
formalism is the following 



S[f;g,ipr, 



2k 2 



d 4 xy^gf(R) + S m [g lll ,,ip m } (1) 



where S m [g, t/j m ] represents the matter action, which de- 
pends on the metric g^ v and the matter fields ip m . For 
notational purposes, we remark that — -^==JJ^-, 
and that the scalar curvature R is defined as the con- 
traction R — g^R^, where R^ v is the Ricci tensor 



R^v — SyX\ v + d\T^ v + r^r^ — r^ p r 



and r^ 7 is the Levi-Civita connection 



(2) 



(3) 



From eq. we obtain the following equations of motion: 



f{R)R llv --f(R)g IJtv -V t N v f(R)+g^Uf{R) = n 2 T^ v 

(4) 

where f'(R) = df/dR. According to eq.(0}, we see that, 
in general, the metric satisfies a system of fourth-order 
partial differential equations. The higher order deriva- 
tives come from the terms V M V„/' and □/'. Only 
when f(R) is a linear function of the scalar curvature, 
f(R) = a + bR, the equations of motion are second-order. 
The trace of eq.(@J is given by 



3D/' + f'R - 2/ = k 2 T 



(5) 



This equation can be interpreted as the equation of 
motion of a self-interacting scalar field, where the self- 
interaction terms are represented by f'R — 2/. This can 
be seen by algebraically inverting the function f'(R) and 
writing R as R — R(f'). In this way, defining 



= /' 

v(4>) = RW - M 



we can write eqs.Q and © as 



-T 



3D0 + 2V{4>)-4>- 



dV 



-9 ni> 



(V M V„<; 



= k T 



(6) 
(7) 



Sm^)(8) 
(9) 



The above equations of motion can also be obtained from 
the following action 



S[g»v, <t>, V>n 



l — J d 4 x^[<j>R(g)- 



2 

w((t>) 



(10) 



(d^d^) - V{cf) 



71 
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which represents a Brans-Dicke-like scalar-tensor theory, 
in the particular case uj = 0. For more details on the 
scalar-tensor representation and a different derivation of 
this result see 181. 



III. COORDINATES AND BOUNDARY 
CONDITIONS 

In order to obtain the metric in the solar system we 
will follow the basic guidelines outlined in chapter 4 of 
Will's book [13 . First we solve eqs.® and © for the 
metric and the scalar field in the cosmic regime, where 
the high degree of homogeneity and isotropy leads to a 
Friedman-Robertson- Walker metric 



With regard to the potential defined for the scalar field, 
see eq.J7J), it is easy to see that dV/d(j> = R |18|. Since, 
the curvature can be expressed as R — Rq + lZ(t,x), 
where TZ(t, x) denotes the local deviation from the back- 
ground cosmic value Rq, it is easy to see that the scalar 
field will not, in general, satisfy the extremum condition 
dV/d<fi = 0. This is to be contrasted with the results 
of the literature regarding the post-Newtonian limit of 
Brans-Dicke-like theories, where it is g enerally assumed 
that the field is near an extremum jlj, [l7| • We thus see 
that for f(R) gravities (or u) = Brans-Dicke-like theo- 
ries) it is necessary to consider all the terms associated 
to the potential. This slight complication, on the other 
hand, has its own advantages, since at the end of the cal- 
culations we may ask Nature about the admissible forms 
that the lagrangian f(R) may have. 



ds z 



g^dx^dx" 



-dt 2 +a{t) 2 dxidx l (11) 



and to a space-independent value of the scalar field, 
= <fi B (t). At smaller scales there appear local devia- 
tions from the cosmic values of the fields. In the solar 
system, for instance, the local deviations will be small, 
thereby allowing us to treat them as a perturbation with 
respect to the background cosmic boundary values g B v 
and (j) B (t). In our computations we will use coordinates 
(i , x 3 ) in which the outer regions of the local system are 
in free fall with respect to the surrounding cosmological 
model. Neglecting the second-order corrections, the 
local and background coordinates are simply related by 
t(to, xq; t, x) w (t — to) and x° (to, xq; t, x) « ao(x — xqY . 
From now on we will omit the bar on the local coordi- 
nates and will denote 0o, 0o the asymptotic boundary 
values of the scalar field at the cosmic time to, i.e., 
0o = S (to) and O = B (i o )- 

For approximately static solutions, corresponding 
to gravitating masses such as the Sun or Earth, to 
lowest-order, we can drop the terms involving time 
derivatives from the equations of motion. In our local 
coordinate system, the metric can be expanded about 
its Minkowskian value as g M „ = ?7 M „ + h^ v . The solu- 
tion for the scalar field can be expressed in the form 
= 0o + i p{ti x )i where 0o = 4>{to) is the asymptotic 
cosmic value, which is a slowly-varying function of the 
cosmic time to, and ip(t,x) represents the local deviation 
from <j)Q, which vanishes far from the local system. We 
want to remark that since 0o and 0o depend on Rq and 
Rq, the metric of the local post-Newtonian system will 
also depend on the background cosmic values Rq and 
Rq. The dependence on these background quantities 
will make the metric change adiabatically in a cosmic 
timescale. This adiabatic evolution could make a theory 
be compatible with the current experimental tests during 
some cosmic era but fail in other periods. We will give 
below some examples to illustrate this effect. In partic- 
ular, we will comment on the Carroll et al. 1 /R model . 



IV. POST-NEWTONIAN METRIC 

As we advanced above, we will expand the equations 
of motion around the background values of the met- 
ric and the scalar field. In particular, we will take 



"'ii.v-i y 



rT ~ hT, = 0o + <p(t,x) 



and V(4>) ~V + ipV ' + <p 2 F " /2 + ... The complete post- 
Newtonian limit needs the different components of the 
metric and the scalar field evaluated to the following or- 
ders 5 oo ~ 0(2) + OU), g 0j ~ 0(3), g tj ~ 0(2) and 
4> - O(2) + 0(4) (see [l2|). The details of the calculations 
and the complete post-Newtonian limit for the theories 
defined in ea. (|10fl can be found in the Appendix. For 
convenience, we will discuss here only the lowest order 
corrections, 500 ~ 0(2), g^ ~ 0(2) and <p ~ 0(2), of 
the case we are interested in, say, ui = 0. The order of 
approximation will be denoted by a superindex. This ap- 
proximation will be enough to place tight constraints on 
the gravity lagrangian. To this order, the metric satisfies 
the following equations 



1 , 
-V 



,(2) 



(2) <y3 



20o + 1 2 0o 20o ) ^ 



= St 



2 

up 
20o 



go 
20o 



20o 



(13) 



where the gauge condition ft£ — fc = dkf^ 2 ' /0o has 
been used. In eliminating the zeroth-order terms in the 
field equation for ip, corresponding to the cosmological 
solution for 0o, the equation for the scalar field to this 
order boils down to 

[V 2 -^]^( £ ,x) = -^ (14) 

where m 2 , is a slowly- varying function of the cosmological 
time given by 



m~ 



00 Vq" 



n 



(15) 
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Note that, despite our notation, there is no a priori re- 
striction on the sign of mz, . The equations of above can 
be easily integrated to give 

^\t,x) 



3 4tt 



d3x 'w^\ Filx ~ x ' l) (16) 



b Air J \x-x'\ 




1 



F(\x-x'\) 



Vo \ \x-x c \ 2 
2</>o J 3 

d 3 x' Pit ' X ' ] 



(17) 



24> 24> 



\x — x'\ 
\x - x c \ 2 



F(\x-x'\ 



+ 



(18) 



where x c is an arbitrary constant vector 1 and the function 
F(\x — x'\) is given by 



-m v \x— x'\ 



if W^p > 



F(\x-x'\) = 



(19) 



cos(m v \x — x'\) if m 2 < 



Note that the term O does not appear in eas. l(TB|l . ((T7|) 
and l|18l) and, therefore, the fact that Ro ma y n °t be 
strictly zero does not affect the Newtonian limit 01 . In 
the post-Newtonian limit it contributes to h 4 Q (see the 
Appendix). In any case, since to all effects </>o is almost 
constant, we can neglect the contributions due to 4>q and 

Since in the solar system the sun represents the main 
contribution to the metric, we can approximate the ex- 
pressions of above far from the sources by 



(2) 

no 

2) 



2G- 



Mc 



60o 



2 

r 



2 lG ^ - ^r 2 
r 60o 



(20) 
(21) 



where M Q = J d 3 x' p SU n(t, x') is the Newtonian mass of 
the sun and the 4>o contribution has been neglected for 
simplicity. We have defined the effective Newton's con- 
stant G as 



G 



1 



F(r) 



87T0O 

and the effective PPN parameter 7 as 

3 - F(r) 



(22) 



7 



F(r) 



1 This vector could be taken, for instance, as representing the cen- 
ter of mass of the system. 



We shall show now that the oscillatory solutions, mr^ < 
— > F(r) = cos(m ¥ ,r), are always unphysical. For this 
case, the inverse-square law gets modified as follows 



M 



1 



cos(m y r) + (m v r) sm(rn y r)^ M Q 



For very light fields, which represent long-range inter- 
actions, the argument of the sinus and the cosinus is 
very small in solar system scales (m^r <C 1). We can 
thus approximate cos(m v r) « 1 and sin(m ¥ ,r) « and 
recover the usual Newtonian limit up to an irrelevant 
redefinition of Newton's constant. However, these 
approximations also lead to 7 ~ 1/2, which is observa- 
tionally unacceptable since 706s ~ 1 19}- If the scalar 
interaction were short- or mid-range, the Newtonian 
limit would get dramatically modified. In fact, the lead- 
ing order term is then oscillating, sin(m„r)M0 /r, and 
is clearly incompatible with observations. We are thus 
led to consider only the damped solutions F(r) — e~ m * r . 

The Yukawa-type correction in the Newtonian poten- 
tial has not been observed over distances that range from 
meters to planetary scales. In addition, since the post- 
Newtonian parameter 7 is observationally very close to 
unity, we see that the effective mass in eas. (|23|) and <|22[) 
must satisfy the constraint m^L 2 3> 1, where L repre- 
sents a typical experimental length scale. Note that when 
u) is not fixed (see the Appendix), there is also the possi- 
bility of having a very light (long-range) field that yields 
almost space-independent values of G and 7. In that 
case, the theory behaves as a Brans-Dicke theory with 7 
given by 



7 = 



l+UJ 



(25) 



and it takes u> > 40000 to satisfy the observational 
constraints pfj| . 

The cosmological constant term (Vo/64>o)r 2 appearing 
in eas. l|20|l and (|21f) also imposes constraints on the par- 
ticular model, since this contribution must be very small 
in order not to modify the gravitational dynamics of lo- 
cal systems ranging from the solar system to clusters of 
galaxies. In the terminology of f(R) gravities, the con- 
straint from the cosmological constant term is 



fo - Rofn 



fo 



Li < 1 



(26) 



where Ll may represent a (relatively large) length scale 
the same order or greater than the solar system. The 
constraint m^Lg 3> 1 associated to the effective mass 



(23) can be reexpressed as 



fo ~ Rafo 
in 
Jo 



4»i 



(27) 



where Lg represents a (relatively short) length scale that 
can range from meters to planetary scales, depending on 
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the particular test used to verify the theory. It is worth 
noting that a generic lagrangian of the form 



f(R) = R + Xh(R) 



(28) 



with A a suitable small parameter, satisfies the two 
constraints of above if h(R), h'(R) or h"(R) are finite 
or vanish as the universe expands. General Relativity, 
which can be seen as the limit A — > 0, saturates those 
constraints. We will consider in the next section some 
examples of theories with the form proposed in eq. (12811 . 
In section IVII we will analyze in detail the implications 
of the constraint of eq. 1(2*71) . 

Before concluding this section, we shall briefly dis- 
cuss some simplifications that we may carry out from 
the above considerations in the complete post-Newtonian 
metric given in the Appendix. First of all, it is worth 
noting that with a tiny Vq we can eliminate part of the 
cosmological constant terms. This fact together with our 
definition for G leads to the PPN parameter (3=1, which 
coincides with the one corresponding to GR. On the other 
hand, a massive field would allow us to neglect the ex- 
ponential terms and the tp( 2 > contributions. Further sim- 
plifications could be achieved from the observational ev- 
idence supporting the constancy of Newton's constant. 
Assuming a massive field, it follows that G/G « —(fio/fio- 
This relation provides a justification to argue that (fio/^o 
and </>o/</>o & r e small, if nonzero. With these simplifica- 
tions we recover the post-Newtonian limit of GR, where 
G/G = 0. We thus see that measurements of a change 
in G with time and of Yukawa-type corrections in the 
inverse-square law could be due to the presence of non- 
linear elements in the gravity lagrangian. 



A. Positive powers of R 

Following the structure of the ansatz proposed in 
ea. (|28|l . we can consider the family of models defined by 
f(R) = R + R n /M 2n ~ 2 , where M represents a very large 
mass scale. We will only consider the cases n > 2. These 
models are characterized by 



V{<j>) 



= M>-l)(^)" = I 2 (n-l)| 



>-r 

n 



(29) 
)(30) 



With eq. (|""T)j) at hand, we can compute the effective mass 
rrv^p that characterizes the post-Newtonian metric. It is 
given by 



Rn 



3(n- 1) 



M2 

3n 



M 2 

Ro 

n-2 



1 



- (n- : 
(n-1)' 



(31) 



where <fio = f'(Ro) and Ro represent the cosmological 
values of <f> and R at the moment t^. The time-time 
component of eq.(["I} can be used to extract some infor- 
mation about the cosmological evolution of R. This will 
help us to understand the adiabatic change in the post- 
Newtonian metric. The expansion factor satisfies the fol- 
lowing equation 



K 2 p - 



R 

M 2 



3rj 



{n — 1) „ „ , . a R 
z an. 



(32) 



V. EXAMPLES 

We shall now illustrate with some simple examples 
how the parameters that define the post-Newtonian 
metric are subject to a slow adiabatic evolution due 
to the cosmic expansion. The aim of this section is to 
point out the relevance of the cosmic boundary values 
of the fields in the description of isolated systems. We 
want to make special emphasis on the fact that the 
gravitational dynamical properties of a local system at 
a given time may not be completely determined by its 
own internal characteristics, but can be affected by the 
state of the universe as a whole at that moment. Only if 
the f(R) lagrangian is linear in R or if the scalar field is 
non-dynamical (Palatini formalism) , the post-Newtonian 
metric is completely determined by the properties of the 
local system. 



Inserting a(t) = doe 7 ' in ea. (|32J) and taking p = 
for simplicity, it follows that at early-times the evolu- 
tion is dominated by the (R/M 2 ) n contribution with 
7 2(n-i) = (M 2 /12)™-7(n - 2). After the early-time in- 
flation predicted by these expansion factors, as the curva- 
ture decays below the scale defined by M 2 , the (R/M 2 ) n 
effect is suppressed and the subsequent evolution is gov- 
erned by GR, with a(t) = a n t s and s = 1/2 during the 
radiation dominated era, and s = 2/3 during the matter 
dominated era. Thus, at all times after the inflationary 
period, we have M 2 /R ^> 1, or equivalently (</> — 1) — > 0. 
This leads to a very large effective mass for the scalar 
field and a tiny cosmological constant term Vo/4>o - > 0. 
In consequence, this family of models yields an accept- 
able weak-field limit. In fact, it seems reasonable to think 
that these theories are compatible with GR in all astro- 
physical applications, since the curvature is expected to 
be much smaller than M 2 in all situations but at the very 
early universe. 
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B. Negative powers of R 

A well-known example of this type is the Carroll et al. 
model H, defined by f(R) = R - fi 4 /R, where [i rep- 
resents a tiny mass scale of order 10 -33 eV. The reason 
for the minus sign in front of /z 4 is intriguing, since this 
definition leads to a negative effective mass 



R 



(33) 



which we have shown to be in conflict with the post- 
Newtonian limit (see ea. l|24|) '). An improved formulation 
of the theory could be obtained by changing the sign in 
front of fi in the definition of f{R). In this way, we can 
easily extend the results of the examples of above to the 
models f(R) = R + fi 2n+2 /R n . A direct consequence of 
the positive sign in front of fj, 2n+2 is the loss of exponen- 
tial solutions for a(t) at late-times, since the relation be- 
tween 7 and fi turns into 7 2 (" +1 ) = -{n + 2)(^ 2 /12) n+1 . 
These models are characterized by 



4> = f = i 



V{cf>) = -fi 2 (n-l) 



n+l 



(34) 



R 



H 2 (n+1) 



(35) 



The effective mass of the scalar field takes the form 



Rn 



3(n 

if! 
3n 



n V A* / 



1 



1+2 
H-l 



(n + 2) 
(n + l)' 



1 



(36) 



We will restrict our discussion to the cases with n > 1. 
The cosmological evolution of these models during the 
radiation dominated era requires a complete solution 
of the model, since a simple power law expansion is 
ill-defined. We will just concentrate on the matter 
dominated era, a(t) — arji 2 / 3 , and beyond, a(t) = dot Sn 
with s n — (2n + l)(n + l)/(n + 2). These solutions 
imply that the curvature decays with the cosmic time as 
R = 6s(2s — l)/t 2 . One can numerically check that the 
transition from the matter dominated era, s = 2/3, to 
its final value s n is smooth (we took K 2 p mo /fi 2 = 3/7). 
During the matter dominated era, n 2 /R — > and <fi m 1, 
eas. (|36|) and (|35|l indicate that m 2 is very large and 
Vo/0o very small. In consequence, these models yield 
a valid post-Newtonian limit. However, as the universe 
expands and the curvature approaches the critical value 
(i? c //x 2 )" +1 = n(n + 2), in which m 2 = 0, the effective 
mass is small and the post-Newtonian limit tends to 
that of a Brans-Dicke theory with oj = 0. At later times, 
m 2 becomes negative and the weak-field approximation 
is ill-defined, as we discussed above. We can, thus, con- 
clude that these theories do not seem a good alternative 



to explain the late-time cosmic speed-up, since they 
have an unacceptable weak-field limit at the present time. 



VI. CONSTRAINED LAGRANGIAN 

A qualitative analysis of the constraint given in 
ca. (l2"7|l can be used to argue that, in general, f(R) 
gravities with terms that become dominant at low 
cosmic curvatures are not viable theories in solar system 
scales and, therefore, cannot represent an acceptable 
mechanism for the cosmic expansion. Roughly speaking, 
ea. (|27|l says that the smaller the term /q , with f({ > 
to guarantee m 2 > 0, the heavier the scalar field 2 . In 
other words, the smaller /q , the shorter the interaction 
range of the field. In the limit /q — > 0, corresponding 
to GR, the scalar interaction is completely suppressed. 
Thus, if the nonlinearity of the gravity lagrangian had 
become dominant in the last few billions of years (at 
low cosmic curvatures), the scalar field interaction 
range would have increased accordingly. In conse- 
quence, gravitating systems such as the solar system, 
globular clusters, galaxies,. .. would have experienced 
(or will experience) observable changes in their grav- 
itational dynamics. Since there is no experimental 
evidence supporting such a change and all currently 
available solar system gravitational experiments are 
compatible with GR, it seems unlikely that the non- 
linear corrections may be dominant at the current epoch. 



Let us now analyze in detail the constraint given in 
ea. (|27|l . That equation can be rewritten as follows 



Rn 



f'(Ro) 
Rof'(Ro) 



- 1 



Li » 1 



(37) 



We are interested in the form of the lagrangian at inter- 
mediate and low cosmic curvatures, i.e., from the matter 
dominated to the vacuum dominated eras. We shall now 
demand that the interaction range of the scalar field re- 
mains as short as today or decreases with time so as to 
avoid dramatic modifications of the gravitational dynam- 
ics in post-Newtonian systems with the cosmic expansion. 
This can be implemented imposing 



Rf"(R) 



> 



1 

PR 



(38) 



as R — > 0, where I 2 -C L| represents a bound to the cur- 
rent interaction range of the scalar field. Thus, ea. (|38|) 
means that the interaction range of the field must de- 
crease or remain short, ~ I 2 , with the expansion of the 



2 Note that </) = /' must be positive in order to have a well-posed 
theory. 
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universe. Manipulating this expression, we obtain 
d\og[f(R)} I 2 



dR 



l + l 2 R 



(39) 



which can be integrated twice to give the following in- 
equality 



f(R) <A + B[R 



PR 2 



(40) 



to that of a cosmological constant and, therefore, do 
not substantially modify the description provided by 
General Relativity with a cosmological constant. 

As a final remark, we want to point out that the 
Starobinsky model f(R) = R + aR 2 [21| besides leading 
to early-time inflation and satisfying the solar system 
observational constraints, also seems compatible with 
CMBR observations E3. 



where B is a positive constant, which can be set to unity 
without loss of generality. Since /' and /" are posi- 
tive, the lagrangian is also bounded from below, i.e., 
f{R) > A. In addition, according to the cosmological 
data, A = — 2A must be of order a cosmological con- 
stant 2 A ~ 10~ 53 m 2 . We thus conclude that the gravity 
lagrangian at intermediate and low scalar curvatures is 
bounded by 



-2A < f(R) < R- 2A 



l 2 R 2 



(41) 



This result shows that a lagrangian with nonlinear terms 
that grow with the cosmic expansion is not compatible 
with the current solar system gravitational tests, such 
as we argued above. Therefore, those theories cannot 
represent a valid mechanism to justify the observed 
cosmic speed-up. 
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APPENDIX A: DETAILED CALCULATIONS 



VII. SUMMARY AND CONCLUSIONS 

In this work we have computed the post-Newtonian 
limit of f(R) gravities in the metric approach using a 
scalar-tensor representation. This representation allows 
to encode the higher-order derivatives of the metric in a 
self-interacting scalar field defined by (f) = df J dR. In this 
manner, the equations of motion turn into a system of 
second-order equations for the metric plus a second-order 
equation for the scalar field. The post-Newtonian metric 
is thus characterized by several quantities related to 
the scalar field, say, the intensity Vq of the potential, 
the length scale ~ m~ l of its interaction, and the 

boundary values </>q and <pQ. Since those magnitudes 
are given in terms of / and its derivatives, we found 
the constraints given in eq. (|2*H|l and eq. ijTH) necessary to 
have agreement between the predictions of these theories 
and the observational data in the solar system. Those 
constraints show that the gravity lagrangian f(R) at 
relatively low curvatures is bounded from above an from 
below according to ea. H41|) . Therefore, f(R) gravities 
with nonlinear terms that grow with the expansion of 
the universe are incompatible with observations and 
cannot represent a valid mechanism to justify the cosmic 
accelerated expansion rate. In the viable models the 
non-linearities represent a short-range scalar interaction, 
whose effect in the late-time cosmic dynamics reduces 



We will take <? M „ « 7 ?M „ + h^, gp» » if - W v 
and <j) — 0o + <p(bi%)- For convenience, we will rewrite 
the equations of motion corresponding to the action of 
eq. (|10|l in the following form 
^2 r 



1 



Ul 



1. 



A2 ^ 



(Al) 



We keep the term with u> because at the same price we 
can compute the post-Newtonian limit of any Brans- 
Dicke-like theory. At the end of the calculations we can 
particularize to the case u> = to obtain the desired 
result. It is also useful to keep the term with oj to check 
that, when the potential terms are neglected, we recover 
the expected limit of Brans-Dicke theories. 

The expansion of the Ricci tensor around the 
Minkowski metric can be written as follows 
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where all the indices are raised and lowered with the 
Minkowski metric. Assuming a perfect fluid, the ele- 
ments on the right hand side of ea. (|All) are given, up 
to the necessary order, by 
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where 
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We can also define the contribution due to the scalar field 
as 



t*, = ^09,0+ iv M V,0+ [D0 + F(0)] (A9) 



Its components are 
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the equations of motion boil down to 
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where simply states the relation hij = Sijh^j}. The 
equation for the scalar field is given by 
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where we have defined 
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